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We show that for appropriate choices of parameters it is possible to achieve photon blockade in 
idealised one, two and three atom systems. We also include realistic parameter ranges for rubidium 
as the atomic species. Our results circumvent the doubts cast by recent discussion in the literature 
(Grangier et al Phys. Rev Lett. 81, 2833 (1998), Imamoglu et al Phys. Rev. Lett. 81, 2836 (1998)) 
on the possibility of photon blockade in multi-atom systems. 



I. INTRODUCTION 

Recent work by Imamoglu et al suggested a promising scheme for observing photon blockade in a highly non- 
linear cavity, where the change in the Kerr nonlinearity due to single photon effects was enough to make the cavity 
non-resonant with modes of more than one photon. As described in their work, such a device could work as a single- 
photon turnstile, similar to that realised recently by Kim et al Q in a semiconductor junction, which might be useful 
for quantum computation and the generation of non-classical light fields. The work by Imamoglu et al was based on 
the use of the adiabatic elimination procedure and later studies HQ showed that the breakdown of the procedure 
in the high dispersion limit leads to prohibitive restrictions on the parameter space where photon blockade could be 
observed in a multi-atom system. Further to this, Werner and Imamoglu || and Rebic et al || suggested that these 
problems could be overcome by using a system with a single atom. The question of observing photon blockade in 
multi-atom systems is of more than just theoretical interest. Possible schemes for observing photon blockade depend 
on directing a low flux atomic beam through a high finesse cavity and this implies that there will be an uncertainty 
in the number of atoms present in the cavity at any given time. We show that the photon blockade, in certain 
circumstances, remains strong even for a fluctuation in the atom number. Recent work |5||jJ has highlighted the 
importance of employing a mutual detuning between the cavity and the semiclassical coupling field to shift the many 
atom degenerate state out of resonance. Our work builds on this idea and presents a detailed atom-cavity dressed 
state calculation showing the parameter regimes which offer the best prospects for the observations of photon blockade 
in one, two and three atom systems. We show that current technology should allow the construction of a system 
exhibiting photon blockade. We suggest that photon blockade can be best realised, not in a MOT as originally 
envisaged by Imamoglu et al [Q, but by sending a low atom current beam through a cavity of the style realised by 
Hood et al B and Miinstermann et al [pi. 

Photon blockade in a cavity can be explained quite simply. An external field drives a cavity that is resonant 
when there are zero or one photon in the cavity and non-resonant for two (or more) photons. This can be achieved 
by introducing a medium into the cavity, exhibiting a large Kerr non-linearity which alters the refractive index 
as a function of the intensity. In general, one would expect to find large non-linearities in systems exhibiting 
electromagnetically induced transparency [|10| . 
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II. DRESSED STATES OF THE ATOM-CAVITY SYSTEM 



Consider a four level atom as depicted in figure 1. The energy levels are labeled in order of increasing energy as |et), 
1 6), |c) and \d) with associated energies of tku a , tuvb, tuvc and ftwd respectively and transition frequencies uj a p — u> a —ujf3 
where a, (3 = a, b, c, d. The \b) — |c) transition is driven by a strong classical coupling field, with frequency w c iass, Rabi 
frequency £1 , detuned from the \b) — \c) transition by an amount S c b — ^ciass — <-^cb- The atoms are in a cavity with 
resonance frequency u cav which is detuned from the \a) — |c) transition by an amount S ca = w C av — w ca , detuned from 
the \b) — \d) transition by an amount A = oj cav — oJdb and not interacting with the \b) — |c) transition. The detunings 
5 c b and 5 ca are set equal to ensure that the \a) — \b) transition is driven by a two photon resonance. We therefore define 
the mutual detuning, 6 = 6 c b = 5 ca . The cavity is driven by an additional classical field with frequency to e = cj cav 
and power, P. We analyse the effect of this field by examining the dressed states of the atom-cavity system. This 
configuration of fields and atomic levels is called the N configuration and has been considered previously 
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The cavity linewidth is r cav . The atom-cavity mode coupling is g ai a 2 = ( w aia 2 /2^£oK;av) l^a 1 a 2 where a\ and a,2 
correspond to atomic levels, fi ai a 2 is the electric dipole moment of the transition, V ca . v is the cavity volume and eo 
is the permittivity of free space. The Hamiltonian for the system with N atoms in the frame rotating at the cavity 
resonance frequency in the rotating-wave approximation is || 

r> N N N 

£ = -if „ £ dl - if d £ d\ d + J2 n [Kb + d L 



N N 

+ 9ac (aa J ca + a ] d 3 ac ) + ^ 9bd (aa db + a^d 3 bd ^j - iT cav a^a. (1) 

J=l 3=1 

where r c = T c + id, Yd = Yd + «A, T a is the decay rate from atomic state \a), a (a'j is the cavity photon annihilation 
(creation) operator and ct^ iQ2 is the atomic operator \a±) (a.2\ acting on atom j. 

We find the dressed states by diagonalizing Ti. Fortunately, H is block diagonal in the bare state basis. The nth 
block can be identified by starting with the state \a, a, . . . ,n), representing all the atoms in state \a) and the cavity 
field in the n photon state \n), and finding the closed set of states coupled to |a, a, . . . , n) by Ti. Diagonalizing this 
block gives the n quanta manifold of dressed states. 

We first consider the case of a single atom in the cavity. The zero quanta manifold consists solely of the state 
\a, 0). The one quantum manifold is spanned by the states \a, 1), \b, 0) and |c, 0) . The corresponding block of H/H 
can be written in matrix form in this basis as 



o o n 

g ac Q -it c 



(2) 



where the superscript on Ti. refers to the number of atoms and the subscript to the number of quanta in the system. 
In order to simplify the expressions for eigenvalues and eigenvectors, we assume r cav = 0. The figures which follow, 
however, have been generated using non-zero values of r cav . Diagonalising the matrix Ti^ /H with T cav = gives 
the dressed state energies 



S + =[-iT c + \l-Tl + A{W+gl c ))/2 
£o = 



and corresponding dressed states \D+), \Dq) and \D-) respectively. In this form, the real part of the eigenstate 
corresponds to the state energy and the imaginary part to the width of the state. These eigenstates form the well 
known Mollow triplet jl2j and are presented in figure 2(a) as a function of the scaled mutual detuning, 5 /CI, with 
r c /fi = 0.1, r cav /0 = 0.01 and g ac /Cl = 1. It is important to express the form of the central dressed state, \Dq) 
which is 

| A,) = ; " \a, 1) - ; 9ac \b, 0) . 
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Photon blockade will occur in this dressed-state picture when the cavity driving field resonantly couples the zero 
to one quantum manifolds, and only weakly couples the one and two quanta manifolds. We can gauge the extent of 
these couplings by treating each transition driven by the cavity driving field as a separate, closed two-state system. 
This approach will break down when multiple states are excited simultaneously. However, in situtations where the 
photon blockade effect occurs, the number of dressed states that are significantly occupied will be minimal and our 
two-state model should give a reasonably accurate picture of the degree of excitation of each transition. 

The effect of the cavity driving field on the cavity-atom system can be treated by including the additional term on 
the right hand side of equation |l| 

h/3{a + a f ) 



where f3 = \J PT CW T 2 / (Ahuj cav ) is the external field-cavity mode coupling strength for a cavity mirror transmittance 
of T. In our two state model the cavity driving field drives transitions between lower and upper states, \L) and \U), 
in the n and n + 1 quanta manifolds respectively. The effective Rabi frequency of the transition is given by 

n e = \(3(L\a\U)\. 

Under these conditions, the steady state population of \U) is given by 

where A e is the detuning of the external cavity driving field from the \L) — \U) transition and Tjj the decay rate of 
state \U), assumed to take population from \U) to \L). Note that we have ignored the decay rate from \L). We 
denote the maximum value of p cxc over all transitions from a given lower state to all possible upper states in the n 
quantum manifold as pcxc- For ideal photon blockade, we require pi]) c ~ 0.5 for the transition from the ground state, 
\L) = \a,a,...,Q) = \Gq) to the maximally coupled one quantum dressed state \Gi). For the case that |Gi) = \D Q ), 
we note that IV will be small, because |D ) contains no proportion of atomic state |c). Thus it is possible to inject a 

(2) 

single quantum of energy into the atom-cavity system for modest values of /?. Ideal blockade also requires that pcxc 
be negligible for transitions between \L) = \Gi) and states \U) of the two quantum manifold. 

We next consider the two atom, one quantum manifold of states. In this case, the basis states are |a, a, 1), \a, b, 0), 
|a, c, 0), |6, a, 0) and |c, a, 0) and the corresponding block of TL/h can be written in matrix form as 



(2) 
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(3) 



(2) 

Diagonalising Tt\ /H yields the eigenvalues 



S+2 = I -*T, 

£+i = (-if, 
£a = 



-f2 + 4(fi2 + 2(&) /2 



-f 2 + 4fi 2 /2 



-f 2 + 4f> 2 /2 



£-2 = ^-*r c - sJ-Tl + 4 (Q2 + 2 5 2 c )J /2 

with associated dressed states l-D+2), l-D+i), \Dq), \D-i) and \D- 2 )- The dressed state energies are plotted in figure 
2(b) for the same conditions as in figure 2(a). There are some important similarities between the spectrum of 
eigenstates for the one atom and two atom cases. In each case there is a state with zero energy, indicating that 
transitions from the zero to the one quantum manifold are possible for a cavity driving field tuned to the cavity 
resonance w cav . The states which are anti-crossing in each manifold are asymptotic to the lines £/Q = and 
£ /n = 6/Q with the point of closest approach being at (5/0 = 0. It is also important to realise that although there 
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are five distinct eigenstates only three of these eigenvalues will couple to the ground state of the atom cavity system, 
i.e. the matrix element (a, 0| a |Djv) is non-zero only for N = 0, ±2. For this reason only the optically active states 
l-D+2), I -Do) and l-D-2) are plotted in figure 2(b). 

We have also solved the analogous three and four atom Hamiltonians and we summarise our results for the eigenstates 
in each case as 



£+2 = -iT c + J-Tl + 4 (ft 2 + Ngl c ) /2 



e+i = (^r c + ^-r 2 + 4ft 2 ) /i 
£0 = 

£-x = (-if c - 7-f2 + 4fi 2 ) /2 



f-2 = -*r c - V -r 2 + 4 (ft 2 + tv 5 2 c ) /2 



(4) 



where TV = 1,2,3,4 indicates the number of atoms in the cavity. The degeneracies of the eigenvalues presented in 
equations |] are interesting to observe, namely the £+2, £0 and £-2 values are all non-degenerate, whilst the £+1 and 
£-1 values are (N — 1) fold degenerate. (The latter values, £+\ and £-\ do not occur for N = 1). The presence of 
the zero eigenvalue, £q, indicates that it is possible to inject one photon into the atom-cavity system with a cavity 
driving field tuned to w C av- 

Now we consider the two quanta manifold for a cavity containing one atom. This manifold is spanned by the states 
|a,2), |6, 1), |c, 1) and \d, 0) and the corresponding block of Tl/h in matrix form is 



— 2fT 

ZjL± cav 





V2g ac 








^r cav 


Cl 






Cl 










gbd 





-if 



In order to observe photon blockade in such a system, it is essential that none of the eigenstates of Ti^ /h are 

resonantly coupled by the cavity driving field to the occupied states of the one quantum manifold. To achieve this, 
(2) (2) 

we require Pex C <C 0.5, the smaller pexc the greater the degree of photon blockade. A plot showing the eigenenergies 
of this system with associated linewidths, T x where x is the dressed state under consideration, is presented in figures 
3(a) and 3(b). In each case A/O = 2, T c /£l — Td/Q — 0.1, r cav /f2 = 0.01, g ac /Q = 1 and in figure 3(a) gbd/Q = 1 
and the energies are plotted as a function of S /Q, whereas in figure 3(b) 6 /CI — and the energies are plotted as a 
function of g^d/Cl. The important features to recognise from these two traces is the shift of the smallest magnitude 
energy state from zero, indicating (as highlighted in || and ||) that photon blockade will indeed be possible in the 
one atom case. 

Similarly, for the two atom case, the basis states of the two quanta manifold are |a, a, 2) , \a, b, 1), \a, c, 1), |&, a, 1), 
|c, a, 1), |a, d, 0), \b, b, 0), \b, c, 0), |c, 6,0), |c, c, 0),and \d, a, 0) where the notation is | atom 1, atom 2, cavity field) . In 
matrix form, the corresponding block of the Hamiltonian is 



o,(2) 
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where T* = T r 
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The energy eigenvalues of Tt 2 /fo are presented in figures 4(a) and 4(b). As in the two atom, single quantum 
case there are optically inactive states, and these have been removed from the figures. Figure 4(a) was generated 
for A/fl = 2, 5/n = 0, T C /VL = Td/Cl = 0.1, r cav /fi = 0.01 and g ac /Q = 1- As can be seen, for this case there 
is a zero eigenvalue so we would expect that there would be resonant coupling from the single quantum manifold to 
the two quantum manifold by a cavity driving field tuned to cij cav and hence, photon blockade would not be observed 
in this case. However in figure 4(b) we illustrate the effect of introducing a small mutual detuning, 6, from state 
\c). The parameters used were A/fl = 2, S/Q = 0.5, T c /fl = T d /fl = 0.1, T cav /fl = 0.01 and g ac /^ = 1- It is 
important to observe the shift of the smallest magnitude eigenvalue from zero. Although this shift is less than for 
the corresponding single atom case, it is still feasible to consider building a cavity to observe this photon blockade. 
This point will be elaborated on in the conclusions p3| . 



III. PHOTON BLOCKADE IN A REALISTIC SYSTEM 



Of critical importance to the experimental observation of photon blockade in atomic systems, is the question of how 
to trap single or small numbers of atoms within a small, high finesse optical cavity. We know of no demonstration of 
continuous trapping, however recent experiments conducted by Hood et al Q and Munstermann et al |9| have shown 
that it is possible to have a cavity with extremely low atom fluxes passing through it. This was achieved in |§| by 
allowing atoms to fall from a leaky MOT into the cavity and in ^ by directing atoms out of a MOT and into the 
cavity. One would expect such atomic ejections to be stochastic in nature and as a consequence, the probability of 
having a certain number of atoms in the cavity should follow Poissonian statistics. In order to build a device which 
uses photon blockade, it will therefore be necessary to ensure that significant photon blockade will be observed over 
as wide a range of atomic numbers as possible, otherwise the photon blockade could be lost and the performance of 
the device degraded. 

A complication in the experimental realisation of photon blockade is that the parameters chosen in the theoretical 
plots shown above were assumed to be independent variables. In general this will not be the case unless great care 
is taken in the preparation of an experiment. Specifically, for transitions in an alkali vapour which might realise 
the N configuration, one would expect the dipole moments of all relevant transitions to be of the same order, with 
the atom-cavity coupling set solely by the cavity volume, a parameter shared by both the \a) — |c) and \b) — \d) 
transitions. One would therefore expect g ac ~ = g. Also, because of the shared cavity, the detuning parameters 
are not independent, so if we assume that u ca — uidb = S u then we find that A = 8 + 5^. These extra considerations 
will be included in the analysis to follow, which concentrates on experimentally realisable effects rather than the 
general theoretical demonstration presented above. We start by investigating photon blockade using the parameters 
Tcav/^ = 0.01, T c /Q, = Td/ft = 0.1, S w /Q, — and f3/Q = 1. For these values pixc ~ 0.5, as is expected due to the 
presence of the strongly absorbing state at the cavity resonance. 

(2) 

In figure 5 (a) we show a pseudo-colour plot of pc* c as a function of 5/D, and g/£l for the one atom two quanta 

(2) (2) 

case, where colour is indicative of the value of Pcxc, blue being and red being 0.5. The graph clearly shows /9 CXC 
decreasing monotonically with g/fl, indicating the effectiveness of photon blockade correspondingly increasing. This 

(2) 

is expected, as a large g will give rise to a highly nonlinear system. Also note that p e x C increases as \5/£l\ increases, 
indicating that the photon blockade is a resonance phenomenon. 

(2) 

In figure 5(b) we show the analogous /?cxc plot for two atoms. In accord with the preliminary results which suggested 
that photon blockade was not possible in the multi-atom system ||f|j|] we find Pexc < 0.5 in the vicinity of 5/Q, = 0, 
implying that photon blockade will not be observable for these choices of parameters. However, by increasing the 
mutual detuning and employing a modest atom-cavity coupling, regions of strong photon blockade are observable, 

/ 2) 

typified by the minimum recorded value on figure 5(b) of pex C = 0.0087. 

The three atom case is shown in figure 5(c) for the same parameter regime as the one and two atom cases. The 

(2) 

form of pexc is similar to that of the two atom case, with the same overall structure, but narrower regions where 
photon blockade should be observable, this is shown by the smaller blue region of figure 5(c) than figure 5(b). Again 

(2) . . . (2) 

very low values of p OX c were obtained, with a minimum recorded value of pc XC = 0.0070. 

To use photon blockade as a tool for new quantum devices it is necessary to identify real systems in which these 
effects may be observed. As an example of what may be achieved with the current state of the art, we take parameters 
from a recent experimental paper |8| and make some minor assumptions about how they might apply to atoms falling 
through a high finesse cavity. It should certainly be possible to achieve a coupling field Rabi frequency of Q = 10 
MHz, external field-cavity coupling strength of (3 = 1MHz and Hood et al achieve g = 120 MHz and r cav = 40 MHz. 
For a system based on transitions in the 87 Rb D2 line |l4| we may assume 5^ = 6600MHz and T c = Td = 17.8 MHz. 
We ignore the Zeeman magnetic sublevels since these form a simple, effective three level A system (for states |a), |&) 
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and |c)) for the situation under consideration [}L5| . Assuming equal rates of radiative decay, these translate into our 
system as g/ft — 12, r c /f2 = T^/Sl = 1.78, r cav /Sl = 4, = 660 and f3/Vl = 0.3. In figure 6(a) we present a plot 

of pixc as a function of g/fl and <5/Sl for one atom in the cavity, with the plots showing pcxc in (b) and (c) for one 
and two atoms respectively. 

There are some important features to note in figures 6. In 6(a) pi^ c increases monotonically with g/Q and, for the 
scale used, is independant of 5/fl. The value of pexl is qualitatively very similar for one, two and three atoms, with 
the values of pix C slightly increasing as the number of atoms increases. As an example, for g/il = 12, pi]} c = 0.3099, 

(2) 

0.3824 and 0.4148 for one, two and three atoms respectively. Traces 6(b) and (c) show pcxc for the one and two atom 

(2) 

cases respectively. The value of pcxc for the one atom case is very small across the entire parameter space, indicating 
that photon blockade should be easy to observe for a single atom. Of interest is the resonance in the vicinity of S = S w 
which will be discussed below. There is an extra resonance in the vicinity of 5 = for the one atom case, although 
this is not present to the same extent in the two and three atom cases. For two atoms, and also for three atoms 

(2) 

(not shown), the overall values (away from the resonances) of pexc appear to increase with the number of atoms and 
are much larger than for the one atom case. Qualitatively, pcxc for two and three atoms are extremely similar, with 

(2) 

maximum values along the line g/fl = 12 of p C x C = 0.2244 and 0.3092 for two and three atoms respectively. These 
observations would appear to agree with the intuitive idea that photon blockade would be more difficult in multi-atom 
systems and that there is a qualitative difference between the single atom case and multi-atom cases. 

(2) 

The results for p cxc in the vicinity of S/il — —660 are shown in figure 7. The significance of this region is that 
we have <5 = — 5 W so that the cavity is now resonant with the \b) — \d) transition. To our knowledge, this situation 
has not been previously explored and the consequences it has for photon blockade are significant. In figure 7 we 
present pixc for one, two and three atoms in plots (a), (b) and (c) respectively. The behaviour of pi^ c shows no 

(2) 

resonance phenomena and is described above, it is only by considering p cxc that the resonance is observed. In 7(a) 

(2) 

we see generally small values for p cxc indicating that photon blockade should be observable. With the exception 
of the 'shelf for g/VL < 1, p cxc has the a similar chevron shape to that observed in the ideal case, shown in figure 
5(a), although with very much smaller values. In figures 7(b) and (c), we observe a roughly triangular region of low 

(2) (2) (2) (1) 

Pexc, superimposed on the background of pcxc noted earlier. Values of pcxc < 0.01 are present where pexc > 0.35 
for both the two and three atom cases. This suggests that the non-linearity of the system is very large about this 
resonance. Further investigations of this resonance will be done with increasing number of atoms to show how robust 
the photon blockade will be as any relaxation of the requirement for low number of atoms will enhance the prospects 
for experimental verification of this effect. 



IV. PHOTON BLOCKADE WITH AN OFF RESONANT CAVITY 



The discussions above have used a system where the cavity is driven resonantly by the external field. This has a 
significant drawback in trying to realise a continuously operating device inasmuch as when there are no atoms in the 
cavity, the cavity will absorb photons from the external field. To overcome this using the on-resonance configuration, 
one must run the experiment in a pulsed mode to ensure that there are no photons in the cavity prior to atoms 
entering the cavity. There is, however, another alternative. This occurs when the cavity driving field is not tuned 
to the cavity resonance, but instead is tuned to a side resonance of the one atom dressed system. In this case there 
will be strong coupling between the external driving field and the cavity when there is one atom in the cavity and 
no coupling when there are zero atoms in the cavity. By studying the form of the dressed state eigenvalues given 
in equation || it is clear that there is a dependence of the eigenvalue with the number of atoms so that it should 
be possible to choose g ac such that the one-quantum manifold is only coupled into when there is only one atom in 
the cavity. These conditions are best met for the case that both 5 and f2 are small compared to g ac to ensure the 
maximum shift in eigenvalues with number of atoms. 



V. CONCLUSIONS 



We have presented detailed calculations which show the parameter space where photon blockade in one, two and 
three atom systems should be observable. By considering transitions on the 87 Rb D2 line, we have suggested that 
photon blockade should be observable using current state of the art technology, such as has been recently demonstrated 
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We have pointed out that strong photon blockade should be possible when the cavity resonance is tuned to the 
1 6) — \d) transition, despite a large mutual detuning, 5. This resonance may relax the requirements on the energy 
levels of the atomic species under consideration. 

We gratefully acknowledge financial support from the EPSRC and useful discussions with Ole Steuernagel (Univer- 
sity of Hertfordshire), Derek Richards (The Open University), Danny Segal and Almut Beige (Imperial College). 



Figure 1: Energy level diagram for the four level N system. The atoms are labelled in order of increasing energy as 
| a), \b), |c), and \d) with energies Hu> a -> hu)b, ftio c and fluid respectively. A strong classical coupling field with frequency 
a-'ciass and Rabi frequency Q is applied to the \b) — |c) transition and detuned from it by an amount 6 = u> c b — w c iass- 
The atoms are placed in a cavity with resonance frequency w cav which is detuned from the \a) — |c) transition by 5 
and from the \b) — \d) transition by A = Udb — w cav . The cavity is driven resonantly by an external classical driving 
field with external field-cavity coupling strength [3 and frequency cj cav . 

Figure 2: Eigenvalues of the one quantum manifold as a function of the scaled mutual detuning <5/f2 for one and 
two atoms in the cavity. In each case the parameters used were g oc /fi = 1, r cav /57 = 0.01 and T c /Q, = 0.1. In figure 
2(a ) the spectrum for a single atom interacting with a single cavity photon traces out the well known Mollow triplet 



JF2|. In figure 2(b) we present the analogous trace for two atoms instead of one where the optically inactive dressed 
states, l-D+i) and |-D_i) have been removed. 

Figure 3: Eigenvalues of the two quanta manifold for a single atom. Parameters used were the same as those in 
figure 2, but with the addition of Td/fl = 0.1 and A/ft = 2. Figure 3(a) has gbd/^t — 1 and the energy eigenvalues 
are plotted as a function of scaled mutual detuning, whilst figure 3(b) has <5 / O = and the energy eigenvalues are 
plotted as a function of gbd/ft- 

Figure 4: Eigenvalues of the two quanta manifold for two atoms in the cavity with the parameters A/ft = 2, 
T c /ft = Td/fl = 0.1, r cav /f2 = 0.01 and g a c/ft = 1, as a function of gbd/ft- In 4(a) the mutual detuning S/ft = and 
no photon blockade is observed, in 4(b) a small mutual detuning of 5 /ft = 0.5 is applied and the photon blockade is 
restored for a cavity driving field tuned to w cav . Note that the optically inactive states have been removed in these 
traces also. 

Figure 5: pcxc plotted as a function of scaled detuning (S/ft) and scaled coupling (g/ft) for one atom (a), two 

(2) 

atoms (b) and three atoms (c). The colour of the plot shows the value of pcx C with blue being zero, red 0.5. The 
parameters chosen for these figures were T cav /ft = 0.01, T c /ft = Td/ft = 0.1, 6 W /Q = and (3/fl = 1. Note that over 
this paramter range, poxc = 0.5. 

Figure 6: Pseudo-colour plots of pi]} c and pcxl- The parameters chosen were those which could be expected in a 
realistic experiment involving 87 Rb atoms. These parameters were T c /0 = Td/fl = 1.78, r cav /fi = 4, S^/ft = 660 
and /3/fi = 0.3. Figure (a) corresponds to pcxc for one atom, whilst (b) and (c) correspond to pexc for one and two 
atoms respectively. Note that the colour scales for each figure vary according to the maximum values of p exc . 

(2) 

Figure 7: Pseudo-colour plots of p?xc in the vicinity of the resonance 5 = — S w for one (a), two (b) and three (c) 
atoms. The parameters used were the same as in figure 6. 
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